##############################
# Table A8
##############################

rm(list=ls())

#sink("~/Dropbox/COVID Peru/08_replication/010_HTE.txt")

library(foreign)
library(Hmisc)
library(ri)
library(stargazer)

##############################
# Prepare data
##############################

# load data
d = read.dta("clean_survey_experiment_covid_peru.dta")

# explore data
names(d)
head(d)
table(d$acuerdo_cuarentena, exclude = NULL)
mean(d$acuerdo_cuarentena_binary)

############################
# HTE: NSE
############################

table(d$nivel2)
d$nivel_3 = as.numeric(d$nivel2)
table(d$nivel_3)

reg1d = lm(d$acuerdo_cuarentena_binary ~ d$t_ind_economia*d$nivel_3)
summary(reg1d)

reg2d = lm(d$acuerdo_cuarentena ~ d$t_ind_economia*d$nivel_3)
summary(reg2d)

############################
# HTE: Informal
############################

reg3d = lm(d$acuerdo_cuarentena_binary ~ d$t_ind_economia*d$informal1)
summary(reg3d)

reg4d = lm(d$acuerdo_cuarentena ~ d$t_ind_economia*d$informal1)
summary(reg4d)

sink()
